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We report on the calculation of the kaon B-parameter using two dynamical flavours of domain wall fermions. 
Our analysis is based on three ensembles of configurations, each consisting of about 5,000 HMC trajectories, with 
a lattice spacing of approximately 1.7 GeV for 16 3 x 32 lattices; dynamical quark masses range from approximately 
the strange quark mass to half of that. Both degenerate and non-degenerate quark masses are used for the kaons. 



1. INTRODUCTION 

The kaon B-paramctcr, Bk, is one of the sim- 
plest quantities that may be calculated on the 
lattice which may be used to directly constrain 
the unitarity triangle. However, its calculation 
on the lattice is greatly complicated when either 
flavour or chiral symmetry is broken by the lat- 
tice action used. When using the domain wall 
fermion lattice action, flavour symmetry is ex- 
act, and the explicit breaking of chiral symme- 
try is small enough to be neglected. Because of 
this it is natural to use domain wall fermions for 
the calculation of Bk, and, indeed, this has been 
done to great success in the quenched approxi- 
mation pp. Here we report on the calculation of 
Bk as part of the first large-scale dynamical (two- 
flavour) domain wall fermion simulations. As well 
as studying the effects of un-quenching, we report 
on the effects of including non-degenerate masses 
(m s ^ m u = ma) in both the lattice calculation 
and the extrapolation to the physical point. 

2. DYNAMICAL ENSEMBLES 

The results of this work are based on three 
ensembles of two-flavour dynamical domain wall 
fermion configurations. Each ensemble utilises a 
lattice of size 16 3 x 32, a domain wall height of 
1.8, a fifth dimension of extent 12, and consists of 
94 configurations with the separation of 50 HMC 
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trajectories between each configuration. The bare 
quark masses of these ensembles are 0.02, 0.03 
and 0.04 respectively, the latter of which roughly 
corresponds to the strange quark mass. Details 
of lattice generation, and further physical results 
can be found in |2I3I4| ; here we will touch on only 
a couple quantities of relevance to the extraction 
of the kaon B-parameter: the residual mass and 
the lattice spacing. All quantities are quoted with 
a jackknife error (although, where relevant, the 
average correlation matrix was used) . 

The residual mass, am rcs , is the additive shift 
in the mass due to the residual chiral symme- 
try breaking of domain wall fermions. We cal- 
culate this from the breaking term in the ward 
identity [5] , extracting our final value of am rcs = 
0.001372(49) by linearly extrapolating to zero dy- 
namical mass. The lattice spacing we use is fixed 
by comparing a dynamical extrapolation of the 
vector meson to the physical point to the exper- 
imental mass of the rho meson. This gives value 
of a^ 1 — 1.691(53)GeV. It is interesting to note, 
however, that should we fix the lattice spacing 
from the Sommer parameter, the result is consis- 
tent 0. 

3. B K 

In the continuum, Bk is defined as 
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where A4 is the time component of the axial cur- 
rent, and 



Oll = S7n(l - 75)rfs7 M (l - 7 5 )d. 
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When using a lattice action which breaks chiral 
symmetry, however, Oll may mix under renor- 
malisation with four wrong chirality operators, 
which we denote O m i Xj i for i = 1, •••,4. Lead- 
ing order chiral perturbation theory predicts that 
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and, unfortunately, that 
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This means that as the chiral limit is approached, 
these wrong chirality operators will dominate the 
calculation. However, for practical lattice calcu- 
lations we do not work in the chiral limit and, 
in fact, we argue that for this calculation the 
mixing with such operators may be entirely ne- 
glected. This may be seen by studying the pat- 
tern of chiral symmetry breaking using the spu- 
rion field technique introduced in 7 ; a simple ap- 
plication of this method leads to the conclusion 
that the wrong chirality operators are suppressed 
by a factor of 0((am res ) 2 ) ~ 10 -6 . A numerical 
study, in the quenched approximation, to bound 
the of size of the contributions of such operators 
has also been made. As reported on elsewhere 
these effects were, indeed, found to be negligible. 

To calculate Bk both the strange quark mass, 
m s , and the average of the up and down quark 
masses, to, must be known. We extract these by 
fitting the quark mass dependence of the pseudo- 
scalar and vector mesons to the predictions of 
NLO partially-quenched chiral perturbation the- 
ory (pq-Xpt) and a linear ansatz, respectively. Re- 
quiring that the pseudo-scalar and vector mesons 
have masses equal to the experimental value for 
the pion and rho in the limit 
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and the pseudo-scalar meson mass be eqaul to the 
kaon mass in the limit 



^dyn — ^valence. 1 — ^ ; ^valence, 2 
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gives a set of equations that we iteratively solve 
to find am = 0.00017(11), am s = 0.0446(29) and 
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Figure 1. (am^) 2 versus bare valance quark mass 
for all ensembles, together with the results of a 
fit to NLO partially-quenched chiral perturbation 
theory. 
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Figure 2. Bare Bk versus timeslice for the 
"idyn = 0.04 ensemble. 

the lattice spacing previously quoted. (Note these 
are bare quark masses; the renormalised quark 
mass is defined as Z m (to + m res ) , where Z m is a 
scheme and scale dependent renormalisation fac- 
tor.) Figure Q] shows the square of the pseudo- 
scalar meson versus (degenerate) valence quark 
mass, together with the results of the fit to NLO 
pq-Xpt • For the calculation of Bk using only de- 
generate valence quarks, we will also make use of 
the strange quark mass fixed by comparing the 
pseudo-scalar mass to the kaon mass in the limit 



^dyn — ^ i ^valence ^s/2 . 
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which leads to a value which is ~ 10% lower than 
the previous one. 

We extract the bare value of Bk from the ra- 
tio given in Eq. [21 using (Coulomb) gauge-fixed 
wall-sources for the kaons at time-slices 4 and 
28, varying the position of the operator. As in 
|5I8) . a combination of propagators with periodic 
and anti-periodic boundary conditions are used 
to avoid contamination due to the periodicity of 
the lattice. Figure [21 shows this ratio for the the 
Wdyn = 0.04 ensemble for various degenerate va- 
lence quark masses. We take our final values from 
an error-weighted average of the points between 
time-slices 14 and 17, but see little change in our 
answer if this range is varied. 

To extrapolate/interpolate to the physical 
point we again use the results of a fit to pq- 
Xpt, restricting the bare quark masses use to be 
< 0.04 in an attempt to stay in the region of va- 
lidity of NLO pq-Xpt- Also, we do not include 
the Wvaicnce = 0.01 data in our analysis due to 
the bad quality of the plateau. Figure [21 shows 
the data for degenerate valence quark masses to- 
gether with the NLO pq-Xpt fit- Completely tak- 
ing the limit given in Eq[Sl leads to a final result, 
when considering only degenerate valence quark 
masses, of 

B d ^ s = 0.509(18) (9) 

while including non-degenerate valence quark 
masses, and taking the limit given in Eq0 gives 

B K = 0.495(18) . (10) 

While these value differ by less than the quoted 
error, as these errors are correlated, this 3% ef- 
fect is statistically well resolved. We quote both 
of these numbers in the MS scheme at 2GeV; 
the needed combination of renormalisation fac- 
tors to convert from the bare value of Bk , Zb k = 
Zll/Z\ J]], is calculated using the NPR tech- 
nique of the Rome-Southampton group. 

4. CONCLUSIONS 

As part of the first large-scale project using dy- 
namical domain wall fcrmions, we have calculated 
the kaon B-parameter in two-flavour QCD. We 
present not only a value using degenerate valence 
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Figure 3. Bare, degenerate, Bk versus valence 
quark mass, together with the results of the fit to 
NLO pq-Xpt- 

quark masses, of B K = 0.509(18) (M5,2GeV), 
but also a value using non-degenerate masses 
(m. ^ m u = m d ), of B K = 0.495(18) 
(MS,2GeV). While this difference is within the 
individual statistical errors, due to correlations 
between the quantities, we find this difference to 
be statistically well resolved. These values are sig- 
nificantly below those found in the quenched ap- 
proximation. We caution that this work is based 
at a single lattice spacing, on a single volume, 
and uses dynamical quark masses which are large 
compared to the physical up and down quark 
masses together with a quenched strange quark. 
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